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1. Introduction 

Currently, K — > nlvi (K/3) decay provides the most precise determination of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element \V US \ through Toe |V Mi | 2 |/ + (0)| 2 , where Tis the decay 
rate and f+{q 2 ) is the form factor defined from the K — > % matrix element of the weak vector 
current 

(n(p')\V^\K(p)) = (p^+p'v) Mf) + (Pn-p' tl ) f-tf), q 2 = ( P -p') 2 - (1-1) 

While the PDG values for the CKM matrix elements [jl|] show a 2a deviation from the CKM 
unitarity 

\V ud \ 2 + \V us \ 2 + \V ub \ 2 = l-8, 5 = 0.0033(15), (1.2) 

recent experiments for T [||, g, |], ||] prefer a slightly larger value of |V„. V |, which is consistent with 
the unitarity (5 = 0). However, in order to make a definite conclusion on this issue, /+(0) has to be 
determined theoretically with an accuracy of about 1%. 

In chiral perturbation theory (ChPT), /+(0) is expanded in powers of meson masses Mk, M„, 
and M n 

f+(0) = 1+/2 + /4 + -", (fn = 0(Ml^)). (1.3) 

Thanks to the Ademollo-Gatto theorem which states that SU (3) breaking effects in /+(0) start 
at 0((m s — m U( j) 2 ), poorly-known low-energy constants of the chiral Lagrangian do not enter the 
ChPT formula for fa [01, and hence it is precisely determined as —0.023. However, this is not the 
case for the higher order corrections /„ (w>4) Therefore, a phenomenological estimate of fa 
based on the quark model [^] has been used in previous analyses of \V US \. 

In order to estimate / + (0) without relying on any phenomenological model, lattice calculations 
have been carried out first in the quenched approximation JTo|], and later in unquenched QCD [11 



12p. In these calculations, however, conventional Wilson- or Kogut-Susskind-type quark actions 
are employed, and hence chiral properties of /+ (0) may be significantly affected by the explicit 
breaking of chiral or flavor symmetry. Since the ChPT formula for fa plays a crucial role in the 
chiral extrapolation of lattice data, it is advantageous to use a quark action, which posses chiral 
symmetry even at finite lattice spacing. In this work, we calculate /+(0) in two-flavor dynamical 
QCD using the DBW2 gauge [13] and the domain- wall quark actions Jl4|] , with which the hadron 



spectrum and the kaon B parameter show good chiral properties [15] 



2. Simulation method 

Our calculations are carried out on a 16 3 x 32 lattice with statistics of 4750 HMC trajectories. 
While we simulate a single value for the lattice spacing cT x =1.69(5) GeV, we expect that the scal- 
ing violation in /+ is not large, since the employed lattice action is (automatically) 0(a) -improved. 
The size of the fifth dimension is set to L5 = 12, which leads to the residual mass of a few MeV. We 
refer to Ref. [ |l5| ] for details on the gauge ensembles used in this study. 

At the moment, we complete our calculation at a single sea quark mass m U( j = 0.02, which is 
roughly half the physical strange quark mass. Three heavier masses m s = 0.03, 0.04 and 0.05 are 
employed for the valence strange quarks. 
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3. Extraction of form factor 

The so-called scalar form factor 



Hq 2 ) = f + (<f) + 



Ml-Ml 



f-(l 2 ) 



(3.1) 



at g max = (Mk — M n ) 2 can be extracted from the double ratio proposed in Ref. [fl6|], 

(M K + M^ 2 



= Cf^/;0,0)Cf (m';0,0) 
l ' ' Cf K (t,t';0,0)C^(t,t';0,0) r,( 



Cj e (f,f';p,p') = £(0 Q (x' ■, OV m (x,04(0,0)) e -'P'("'-") e -^, = * or*), (3.3) 



4M i? M 7t 



|/ (<4Jn (3.2) 



max) 



where represents the interpolating operator for pion (kaon). As shown in Fig. [j], fo{q\ 

is determined with an accuracy of < 0.1 %, since various uncertainties of the three-point function, 
such as the statistical fluctuation, are canceled at least partially in the ratio. 
To study the q 1 dependence of the form factor, we calculate 



F (p,p') 



ttA f 1 + l^lffl « <»") ) ■ « (»") = /- <«">//♦ ' 1 



from a ratio 



Cf*(t,t'; P ,p')C K (t;Q)C*(t'-t;Q) 



E K (p)+E n (p') 
M K + M n 



F(p,p'), 



(3.4) 



(3.5) 



Cf n (t , t' ; 0, 0) C K (t ; p) C n (t' - 1 ; p') f,(*'-0^» 

where C n ( K '(t;p) represents the pion (kaon) propagator with the spatial momentum p. Figure |j] 
shows a plot of F(p,p'). The accuracy of F{p,p') is typically 5-10% for the spatial momenta 

IpUp'I < V2. 

In order to convert F(p,p') to fo(q 2 ), %(q 2 ) is estimated from the double ratio proposed in 
Ref. |0] 



^ 2 ) 



(*=1,2,3), 



Cf^t,t';P,P')C§ K (t,t';P,P') 
(E K (p)+E K (p'))(p+p') k + (E K (p)+En(p'))(p+p') k R k 



(E K (p)+E K (p>)) (p-p% - (^(p)-^(p')) {p+p') k R k 



(3.6) 
(3.7) 



a m s =0.03 
■ in=0.04 
A m =0.05 




1 { i iMiinnm ■ ■ a • 



lpl=l, lp'l=0 
#m =0.04 



Figure 1: Plots of /o(<?m ax ) fig ure ) an d F(p,p') with |p| = 1 and |p'| =0 (right figure) as a function of t . 
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We observe that ^(q 2 ) has a mild dependence on the valence quark mass and its magnitude is 
typically ~ -0.01 with 50-100% error. 



4. Interpolation to q 2 = 

Here, we test two methods to determine /+(0) (=/o(0)). 



method- 1: As in Ref.[10], we first calculate fo(q 2 ) at each simulated q 2 from fo{q 
F(p,p) and £ (q 2 ), and then interpolate fo(q 2 ) to q 2 =0. 



1 ) 

max/ 



method-2: As in Ref.[[Tl|], F(p,p') (%(q 2 )) is interpolated (extrapolated) to q 2 = 0, and then 
/o(0) is calculated at q 2 =0. The interpolation of F (p,p') is carried out using data with fixed 
|p| (or |p'|) so that we can unambiguously identify |p'| (|p|) corresponding to q 2 = 0, which 
is needed to convert F(p,p')\ q 2 =0 to /o(0). We repeat this analysis for two data sets with 
p = and p' = 0, and take the average of results for /o(0). 



We test quadratic and polar fits 



2 ^ = ff{0)-{\+ Cl q 2 +c 2 q A ) 

2^ / " 



1 - c\ q 2 



(4.1) 
(4.2) 



for the interpolation of fo in method- 1 and F for method-2, while only the quadratic fit is tested 
for We observe that four possible ways (quadratic or polar fit for method- 1 or 2) give mutually 
consistent results with an accuracy of < 1%. This is because we have very accurate data of /o(<7max) 
near q 2 = 0, and hence the uncertainty due to the choice of the interpolation method is not large. 

However, we observe that method-2 leads to a slightly smaller error of /o(0) than method- 1, 
since the kinematical factor (Ek{p)— E J[ (p'))/(Ek(p)+E J[ (p')) in Eq. ( fDt ) is not large at q 2 = 0, 
and hence the uncertainty of E, (0) has small influence to /o(0) in method-2. Therefore, we employ 
/+ (0) obtained from method-2 with the quadratic fit in the following. 
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Figure 2: Interpolation of fo(q 2 ) (left figure) and F (p,p') with p' = (right figure) to q 2 — 0. Both figures 
show data at m s =0.04. 
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5. Chiral extrapolation 

From the Ademollo-Gatto theorem, the higher order correction 

00 

A/ = /+(0)-(l+/ 2 ) = £/2* (5.1) 

k=2 

is proportional to (m s —m u d) 2 . Therefore, as in Ref.JlO|], it is convenient to consider a ratio 



for the chiral extrapolation of /+(0). In this analysis, fi_ at the simulated quark masses is calculated 
by using the ChPT formula in unquenched QCD [{7|] 1 . 

Our results for R&f are plotted in Fig. ||[ Since the single sea quark mass is simulated so 
far, we extrapolate R&f only in terms of the valence quark mass by naively assuming that its sea 
quark mass dependence is not large below the simulated sea quark mass. From a simple linear 
fit shown in Fig. |3[ we obtain /+(0) =0.955(12), which is consistent with both of the previous 
phenomenological estimate 0.961(8) [§] and the quenched result 0.960(9) [JTot] - We note that recent 



unquenched calculations in Refs.Jil|, |T2| ] also obtained similar values for /+(0). 
6. Conclusions 

We have calculated /+(0) in two-flavor QCD using the domain-wall quarks. Our preliminary 
result is consistent with the previous phenomenological estimate with an accuracy of 1% level. 
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Figure 3: Quark mass dependence of R&f. Open circle and square represent values corresponding to the 
phenomenological estimate Q and the quenched result [ 10 1. 



'Very recently, fi is calculated in partially quenched ChPT (PQChPT) J14]. We confirm that /+(0) obtained from 
the chiral extrapolation with the PQChPT formula for fi is consistent with the result presented in the text within the 
statistical errors. 
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Our estimate combined with recent experimental results of T leads to |V MS | which is consistent 
with the CKM unitarity. For instance, by using T from the E865 Collaboration [gj, we obtain 
\V m \ =0.229(4) which leads to 

\Vud\ 2 + \Vus\ 2 + \Vub\ 2 = l-8, 8 = -0.001(2). (6.1) 

We note that, however, our result has an additional uncertainty arising from the fact that we have 
not taken the limit of the physical sea quark mass. To remove this uncertainty, our simulations at 
two different sea quark masses are in progress. 

The nice consistency between our and previous quenched estimates of /+ (0) may suggest that 
the systematic error due to the quenched approximation for strange quarks is not large. This point, 
however, has to be confirmed by extending our calculation to three-flavor QCD. 

We thank RIKEN, BNL and the U.S. Department of Energy for providing the facilities essen- 
tial for this work. The work of TK is supported in part by the Grant-in- Aid of the Japanese Ministry 
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